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ABSTRACT 



In this paper we compare Adomian decomposition method (ADM) and Shifted Chebyshev polynomials method in 
order to obtain an approximate solution of nonlinear fractional integro-differential equations of Volterra and Fredholm 
integro-differential equations .We present some examples to find out accuracy of each method. 
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1. INTRODUCTION 

In recent years, there has been continuously renewed interest in integro-differential equations. Many mathematical 
models of physical phenomena produce integro-differential equations as fluid dynamics, biological models, and chemical 
kinetics, ([1],[3],[5]). Nevertheless, the development of the theory of derivatives and integral is due to Euler, Liouville and 
Abel (1823). However during the last ten years fractional calculus starts to attract much more attention of physicists and 
mathematicians, ([7], [11], [14]). Some scientists and researchers interested in searching for method to approximate 
numerical for getting solution of the integro-differential equations of order fractional such as, Adomian decomposition 
method has been widely used by many researchers to solve the problems in applied sciences (Adomian 1944; Adomian 
1989; Kaya and El-sayed 2003). Decomposition method provides an analytical approximation to linear and nonlinear 
problems. 

In this method the solution is considered as the sum of an infinite series, rapidly converging to an accurate 
solution Shifted Chebyshev polynomials is applied for solving fractional integro-differential of equations the following 
form: 



for x; t 2 [0; 1], X is a numerical, where the function g(x), k(x; t) are known and y(x) is the unknown function, D a 
is Caputo fractional derivative and a is a parameter describing the order of the fractional derivative and F(y(x)) is a 
nonlinear continuous function. 




(i.i) 



And 




a > 0. X > 0 . 



(2.1) 
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2. BASIC DEFINITIONS 

In this section we present some basic definitions and properties of the fractional calculus theory, which are used in 
this paper. 

2.1 Definition 

The Rieman-Liouvill fractional of a >0 is defined as, [8]. 



I a //(•'•) = / ( x ~ 0 a_1 2/(*) « > 0, x > 0 

It has the following properties: 

I°y(x) = y(x). 

Where 10 = I (Identity operator) 
2.2 Definition 

The Caputo Definition of fractional derivative operator is given by: 



D a y{x) = 



J m - a y m (x), m-l< a < m . m € N 

D m y(x) 



It has the following properties: 

D a I a y(x) = y(x\ 



■m — 1 h. 



(2.1) 



(2.2) 



(2.3) 



(2.4) 



(2.5) 



I a D a y(x) = y(x) - !/ W (0+) — , x>0. 

k=0 (2.6) 

r(7 + l - a) 

Forx>0; a>0; y>-l 

Caputo fractional differentiation is a linear operation, similar to inter order differentiation. 

D a [Xy(x) + ng(x)] = XD a y(x) + [iD a g(x), where A and /i are constants,( [10] , [15] 



(2.7) 
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3. NUMERICAL SCHEME 

In this section Adomian decomposition method and Shifted Chebyshev polynomials are applied for solving 
nonlinear fractional integro-differential equations. 

Adomain Decomposition Method 

Consider equations (1.1) and (1.2) where D a is the operator defined by (2.4) and (2.1) operating with T on both 
sides of equations (1.1) and (1.2) with obtain: 



D a [Xy(x) + w (x)] = \D a y(x) + ^D a g(x). 

And 



m—1 



71=0 

And the nonlinear function F is decomposed as: 

oo 

' v 



f = J> 



n=0 

Where, A n are the adomian polynomials given by: 



1 d n 



(3.1) 



//(.'•)= ^// (A-) (0 + )^--P(//(.r)-A / k(x r t)F{y(t))d^ 

k=0 J ° (3.2) 

Adomain method defines the solution y(x) by the series, ([2], [12], [13]) 



(3.3) 



(3.4) 



1=0 (3.5) 
The components y 0 , yi, y2 are determined recursively by: 

m— 1 ^ 

2/o(.<-)= ^2/ W (0 + )^r + / a (yW), 

fc =° (3.6) 

y k+1 U)=I a (g(x))+I a { I* kU:.t)A k (t)dt). 

^° (3.7) 
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y k+1 (x) = r(a(x)) + /»( / k(xJ)A k (t) dt). 

^° (3.8) 

Having defined the components y 0 , y l9 y 2 . . . the solution y in series form defined by (3.3) follows immediately. 

Shifted Chebyshev polynomials method 

Shifted chebyshev polynomials is applied to study the numerical solution of nonlinear fractional integro- 
differenetial equations. 

This method is based on approximating the unknown function y(x) as: 

n 

i=Q (3.9) 
a i? i=l,2.... are constants. 

Where Ti (x) is the shifted chebyshev polynomials of first kind which is defined in terms of the Chebyshev 
polynomials T n (x) by the following relation, ([6], [9]). 



t:(.v) = t:(2.l- - 1) , 

And the following recurrence formulae: 

T*(x) = 2(2,: - mt-M ~ , n = 2,3... , 

With the initial conditions 

T 0 *(.r) = l , T?(x) = 2x-l, 

And the collocations points 

% = + cos { — )] * 3 — ^ 1, ... n . 
2 n 

Substituting from (3.9) into (1.1) and (1.2) we obtain 



(3.10) 



(3.11) 



(3.12) 



(3.13) 



D a (jr i a i T?(x)) = g(x) + \ f k(.r.t) M 

i=0 Jo i=0 



(3.14) 



and 
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i=0 J ° i=0 (3.15) 

Equation derived from (3.14) and (3.15) can be written as 
n rX n 



Similarly 

n «i n 

? = 0 *° 2=0 

Substituting from (3.13) into (3.16) and (3.17) we have: 

j2«iD a T*(.c j ) = g( Xj ) + A f ' k(.r j .t)CTa i T;(t))<lt . 
i=o Jo i=o 

And 

p -»1 n 

1=0 ^° z=0 



(3.16) 



(3.17) 



(3.18) 



(3.19) 



From (3.18) or (3.19) we obtain system of nonlinear algebraic equations, solving this system we obtain the values 
of the constant ao, a b a 2 substituting from these constants into (3.9), we obtain: 

y(x) = a 0 TS(x) + aiT?(x) + a 2 T£(x) + ... . (32Q) 

4. NUMERICAL EXAMPLES 

In this section we present some numerical example of nonlinear fractional integro-differential equations adomian 
decomposition method and shifted Chebyshev polynomials and compare the results. 

Example 4.1 

Consider the following nonlinear fractional integro-differential equation: 

D°' 7B y(x) = g(x) + f xt(y{x)f dt, 

J ° (2.1) 

And 
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, , 64 j; 9 / 4 V2r(3/4) x 

g(x) = — -, 

15 ~ 8 



With the initial condition y(0)= 0 and the exact solution y(x)= x 



The Solution According to (ADM) 

y 0 (x) = y (0) + I 3 ' i (g(x)), 



k=0 



m ( x ) = x 3 - 0.7771894662 .r 7 / 4 . 

//!(•'•) = / 3/4 ( [ xt(A Q {t))dt, 
Jo 

yi Or) = 0.06408417608 * 7/4 , 



y n +i(*) = 7 : 



^([\A n (t))dt. 
Jo 



(4.2) 



(4.3) 



(4.4) 



(4.5) 



Then y ^ - yo( x ) + yi( x ) + - < 

The Solution According to Shifted Chebyshev Polynomials 

y{x) .r 3 - 0.01363477054 x 7 / 4 , 

When part of the truncated the summation in (3.9), taken n=3 and substituting in (4.1) we obtain 
Impact Factor(JCC): 4.3912 NAAS Rating : 2.72 
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15 7T 8 J 0 



(4.6) 



By using the collocation point which is defined in (3.13) we obtain system of nonlinear algebraic equations: 



D 3 / 4 (^a n r 7 :(,-)) = - — V 7r w - g + / *t(5>„r*(*)) 2 <ft, 

^° 77=0 



77=0 



(4.7) 



Solving this system we obtain the values of the constant a 0 = 5/16, ai=3/16, a 3 =l/32, substituting from those 
constant into (3.20). 

V ri*/*T*l S 64 -r; /4 ^2r(3/4) , f » 

2^«» £> ' ^(./-j) = — — + / *j«(2^ ,h - 

n=0 *' u ?i=0 

We obtain the approximate solution of equation (4.1) which is the same as the exact solution. Table 1 and figure 1 
shows the comparison between shifted chebyshev polynomials method and (ADM) method. 



15 



3 



»(*) = ^ + g(2x - 1) + -(8./^ - 8x + 1) + ^(32x 3 - 48x a + 18.r - 1) = x* . 




approximate of (ADM) 

approximate of Shifted Chebyshev Polynomials 



Figure 1: Numerical Results of Example 4.1. 
Table 1: Indicate the Amount of Error in Example 4.1. 
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X 


exact = sh 1 1 f < d c}}( ' l)i/slic v polynomials 


Approximant by^ADAI ) 


Error of (ADM) 


0.1 


0.001 


0.00075755356829 


0.0002424643171 


0.2 


0.008 


0.007184450500 


0.000815549500 


U.o 


n no i 
U.Uz / 


U.Uzoo41yUz04 




0.4 


0.064 


0.06125682940 


0.00274317060 


0.5 


0.125 


0.1209463585 


0.0040536415 


0.6 


0.216 


0.2104228475 


0.0055771525 


0.7 


0.343 


0.3356958547 


0.0073041453 


0.8 


0.512 


0.5027731107 


0.0092268893 


0.0 


0.729 


0.7176610661 


0.0113380339 


1 


1 


0.9863652295 


0.0136347705 



Example 4.2 

Consider the following nonlinear fractional integro-differential equation, [4] : 



8x 15 256x 45 



With the initial condition y(0)= 0 and the exact solution is y(x)= x . 
The Solution According to (ADM) 



D^iAx) = 7 



S,3/2 



256.r 9 / 2 



+ i 0 - 5 (y(tyy 



3r(l/2) 315r(l/2) 

Using equation (2.1) in (4.9) we obtain: 



, /2 , v 8.,- 3 / 2 256,- 9 / 2 
D 1/2 y(.r) = ^tt^t - ^tttt^t - 



3r(i/2) 3i5r(i/2) r(i/2) 



mi (l 



fr 1/2 (!/(f)) 2 <« 



m— 1 



//(>) = E* (fe Ho + % + / 1/2 (: 



8.c 3 / 2 256.r 9 / 2 



k=0 



3r(l/2) 315r(l/2) 



- 



/1/2( T(1/2)1 



(4.8) 



(4.9) 



(4.10) 



f)- 1 / 2 (. (/ (f)) 2 *). 



"il 1 8 .r 3 / 2 256.r°/ 2 

.Vo(.'-)= E-'/' < U )7T + / / 



3r(l/2) 315r(l/2)" 



(4.11) 



(4.12) 



yo(x) = -0.2x s +x 2 , 

1 P'JC 

>nu) = /1/2 (f(i7^ / o (•'• - tr 1/2 (Ao(t))rit) . 



(4.13) 
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y x {x) = -0.05a; 8 + 0.2.r 5 + 0.0036363G3G35.r n . 
V2(x) = ^ 1/2 (^T f (■'• - T 1/2 (Ai(0) dt) . 

U1/ - ) - y ° (4.14) 
m (x) = 0.05./ 8 + 0.00008556149734., 17 - 0.01636363636.C 11 + 0.001 94805 1948.r 14 . 

y« + i(x) = Il/2 (f^j£( x ~ t)- 1/2 (A n (t))dt) . 



(4.15) 



Then 



y(x) = yoU) + yi(.r) + y 2 U) + ... , 

y{x) * -0.00008556149734.r 17 + 0.001948051948a- 14 - 0.01272727272X 11 + x 2 , 

Is the approximate solution 
The Solution According to Shifted Chebyshev Polynomials 

Using equation (2.1) in the third reduction of right side in (4.9) we obtain: 

D 1 " tj ) = mm " 3i 5 r(i/2) + ru72) L " - t] m) *■ 

When part of the truncated the summation in (3.9), take n=2 and substituting in (4.1 1) we obtain 

2 



(4.19) 



(4.10) 



(4.16) 



By using the collocation point which is defined in (3.13) we obtain system of nonlinear algebraic equations 

2 Bx m 256 r 9/2 1 f*J 2 



Solving this system we obtain the values of the constant a 0 =8/3, ai=l/2, a 2 =l/8, substituting from those constant 
into (3.20). 
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.'/<•'•) = | + §{2* - 1) + |(8* 2 - 8x + 1) = ,- 2 . 

We obtain the approximate solution of equation (4.9) which is the same as the exact solution. Table 2 and figure 2 
shows the comparison between shifted chebyshev polynomials method and (ADM) method. 




approximate of (ADM) 

approximate of Shitted Chebyshev Polynomials 



Figure 2: Numerical Results of Example 4.2. 
Table 2: Indicate the Amount of Error in Example 4.2 



y 


exact. = shifted chebyshev polynomials 


Approj-rrt m i if by [ A DAI] 


Error of (ADM) 


OJ 


E.UU 


0.01 


0 


0,2 


0,04 


0,3999999974 


2.6 x 10" 10 


0,3 


0,09 


0.08999997754 


2,245 x 10" e 


0.1 


0.16 


0.151)991)4711 


x io- 7 


0,5 


0,25 


0,2499939037 


i\M\m w.m 


0.G 


0,3G 


0,3599553379 


0,0000440021 


0.7 


0,49 


0.4^)70 13633 


(Ul002:Mi4ft7 


0.8 


0.64 


0,6389904851 


0. 0010095149 


{).[) 


0.H1 




0JX)35(i25704 


J 


1 


0.9S&1352177 


0.010SG17823 



Example 4.3 

Consider the following nonlinear fractional integro-differential equation: 



4 j 



0.75 



3r(0.75) 

With the initial condition y(0)= 0 and the exact solution y(x)= x 
Impact Factor(JCC): 4.3912 
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(4.18) 
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The solution according to (ADM): 

m— 1 a 3/4 /*1 

//<.'•) = g^) ( o+)^ + /V4 ( |_^__| )+/ i/4 ( jf 

m — 1 /, /J 3/4 

k=0 v ' 7 

«,„(.,•) = -0.2206525301.r 5 / 4 + 0.9999999999.1- . 

.</!(.'•) =I 1 ' A ([\-t{A 0 (t))dt). 
Jo 

yi(.r) = 0.13S55475G7.r 5/4 . 
Jo 

y 2 {x) = 0.0455S54B20aa; 5/ *. 

2/ 3 (s) = I l / 4 ( / xt(A 2 (t))dt), 
Jo 

y 3 (x) = 0.01874349104.r 5/4 . 

<ln + lU)=I 1 / 4 ( J' .ct(A n (t))dt). 

Jo 

T/nn ! j( J -)^y 0 (r)- f j 1 (r)^y 2 (x)+y 3 ir) + ... . 

y(x) ^ 0.9999999999.r - 0.01779881946a 5 / 4 , is the approximate solution. 
The Solution According to Shifted Chebyshev Polynomials 



a ..3/4 r rl 



(4.19) 



When part of the truncated the summation in (3.9), taken n=2 and substituting in (4.19) we obtain 

2 4 r 3 / 4 r f 1 2 

£ 1/4 (X>^*(.r))= 37^777; -T + / 

' ° (4.26) 
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By using the collocation point which is defined in (3.13) we obtain system of nonlinear algebraic equations, 



2 I r 3 / 4 ,. rl 2 

n=0 "W 4 ! 4 JO ?1 = 0 

(4.27) 

Solving this system we obtain the values of the constant a 0 = 1/2, ai= 1/2, a 2 = 0, substituting from those constant 
into (3.20). 




approximate of (ADM) 

approximate of Shifted Chebyshev Polynomials 



Figure 3: Numerical Results of Example 4.3. 

Table 3 and figure 3 shows the comparison between shifted chebyshec polynomials method and (ADM) method. 



Table 3: Indicate the Amount of Error in Example 4:3 



./' 


exact = shifted chebyshev polynomials 


Appi'o.rimanl by(ADM ) 


Error of (ADM) 


0.1 


0.1 


0.9899909882 


0.00100090118 


0.2 


0.2 


0.1976194424 


0.0023805576 


0.3 


0.3 


0.2960482199 


0.0039517801 


0.4 


0.4 


0.3943380480 


0.0056619520 


0.5 


0.5 


0.4925165183 


0.0074834817 


0.6 


0.0 


0.5906010300 


0.0093989700 


0.7 


0.7 


0.6886037032 


0.113962968 


0.8 


o.s 


0.7865335326 


0.0134664674 


0.9 


0.9 


0.8843974948 


0.0156025052 


1 


1 


0.9822011804 


0.0177988196 
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CONCLUSIONS 

From solving examples we find that Shifted Chebyshev polynomials method is better than Adomian 
decomposition method, the results obtain by using Maple 16. 
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